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We use a simple real-space renormalization group approach to investigate the critical behavior of 
the quantum Ashkin-Teller model, a one-dimensional quantum spin chain possessing a line of criti¬ 
cality along which critical exponents vary continuously. This approach, which is based on exploiting 
the on-site symmetry of the model, has been shown to be surprisingly accurate for predicting some 
aspects of the critical behavior of the quantum transverse-field Ising model. Our investigation ex¬ 
plores this approach in more generality, in a model where the critical behavior has a richer structure 
but which reduces to the simpler Ising case at a special point. We demonstrate that the correlation 
length critical exponent as predicted from this real-space renormalization group approach is in broad 
agreement with the corresponding results from conformal field theory along the line of criticality. 
Near the Ising special point, the error in the estimated critical exponent from this simple method is 
comparable to that of numerically-intensive simulations based on much more sophisticated methods, 
although the accuracy decreases away from the decoupled Ising model point. 


I. INTRODUCTION 

Real-space renormalization group (RG) methods have 
a long and successful history in the study of critical be¬ 
havior in quantum many-body systems in El- Research 
into such methods remains very active, as new innova¬ 
tions seek a delicate balance: to be tractable, the expo¬ 
nential size of the state space for large quantum many- 
body systems requires an extreme truncation of degrees 
of freedom, yet discarding too much can lead to inac¬ 
curate predictions. Real-space RG methods beyond the 
original blocking approach of Kadanoff [3] include the 
density matrix renormalization group [4], projected en¬ 
tangled pair states [5], the multiscale entanglement renor¬ 
malization ansatz (MERA) [6], and tensor network renor¬ 
malization 7], among others. 

In contrast to these sophisticated techniques, it has 
recently been shown that a very simple real-space RG 
map can be used to accurately predict some of the 
critical behavior of the quantum transverse-field Ising 
model [8, 9, Rather than requiring intensive numeri¬ 
cal calculations, this method can be solved analytically 
and gives closed-form expressions for RG fixed points and 
critical exponents. This simple method exactly predicts 
the correlation-length critical exponent v = 1 for the 
one-dinrensional model, and is surprisingly accurate for 
higher-dimensional models on a variety of lattices M- 
The method is shown to yield accurate results for the 
Potts model as well [9,. It is unclear, however, whether 
this success is particular to the Ising model (and variants, 
like the Potts model) or can be applied to a wider range 
of models and more general studies of criticality. 

In this paper, we formalize and generalize this simple 
real-space RG technique to preserve on-site symmetries 
under blocking; we refer to this method as the symmetry- 
respecting real-space renormalization group (SRS RG). 
We then apply the SRS RG to the quantum Ashkin-Teller 
(AT) model [T0| (TTj. This one-dimensional quantum spin 
lattice model has several properties that make it desir¬ 


able for an investigation into the general efficacy of the 
SRS RG map. First, it includes many features of the 
Ising model, such as the on-site symmetry that has been 
argued to be a key component of the success of this sim¬ 
ple technique Hi- In fact, for a particular choice of 
Hamiltonian parameters, the AT model reduces to two 
uncoupled Ising models. The AT model has a much 
richer structure, however. Most notably, it possesses a 
line of criticality, along which the correlation length crit¬ 
ical exponent is expected to vary continuously. There¬ 
fore, calculating this critical exponent using the SRS RG 
along this line is a useful test of the method’s perfor¬ 
mance. We note that the AT model along this line of 
criticality is equivalent to the exactly-solvable six-vertex 
model [12], and is expected to be described by a con¬ 
formal held theory (CFT) with central charge c = 1, 
specifically the so-called c = 1 orbifold boson CFT. We 
can compare the calculated critical exponents from SRS 
RG with those obtained from the six-vertex model [12] as 
well as the known behavior of the CFT. A recent study of 
the critical behavior of the AT model and a detailed com¬ 
parison with the c = 1 orbifold boson CFT predictions 
was performed using MERA [l T5] (a much more sophisti¬ 
cated and numerically-intensive real-space RG method), 
and this provides a useful benchmark. 

Our analysis shows that this symmetry-respecting real- 
space RG scheme can quantitatively capture aspects of 
the critical behaviour for the AT model, and performs 
reasonably compared with more sophisticated numerical 
methods especially given its simplicity. By focussing on 
preserving the symmetry at each blocking, the SRS RG 
captures many of the important features of the RG flow 
and behavior along the critical line. Although this tech¬ 
nique may be of limited use as an accurate numerical sim¬ 
ulation method, its simplicity and analytical form may 
allow this method to play a role analogous to mean-field 
theory for critical RG flows. 

Specifically, when applied to the AT model, we find 
that the SRS RG identifies fixed points on the critical 
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line. Although the scaling field associated with flow along 
this line is predicted from the CFT to be marginal, we 
find that the SRS RG scheme does not exactly repro¬ 
duce this marginal behavior and instead yields two RG 
fixed points on the critical line. Nonetheless, the RG 
flow can be used to calculate the correlation length criti¬ 
cal exponent along the entire line of criticality, and gives 
a simple closed form expression for this exponent that 
is in broad agreement with the CFT behavior. The cal¬ 
culated exponent is exact at the Ising decoupling point 
(the point in parameter space where the AT model re¬ 
duces to two uncoupled Ising models) and the relative 
error in this calculated exponent is comparable to the 
numerically-intensive MERA simulations performed in 
Ref. [13]. Specifically, the SRS RG gives a relative error 
in the prediction of < 10% in the range A € [—0.4,0.9], 
and much small error for |A| -C 1, compared with < 5% 
error in Ref. [13]. The accuracy decreases significantly 
towards the endpoints of the critical line. While the sim¬ 
plest SRS RG is obtained using the smallest blocking 
size of 2 to 1, we also explore alternate larger blocking 
schemes, which yield similar results. 


Our paper is structured as follows. In Sec. [IT] we intro¬ 
duce the SRS RG scheme in terms of a tree tensor net¬ 
work, focussing on the construction of the basic isometry 
defining the RG map through properties of the Hamil¬ 
tonian and its on-site symmetries. In Sec. Ill we apply 
the SRS RG to the quantum Ashkin-Teller model, con¬ 
structing this RG map and using it to determine the fixed 
points and critical behavior of this model along its line 
of criticality. We compare the results to the predictions 
obtained from the orbifold boson CFT, and also explore 
alternate blockings. We conclude and discuss these find¬ 
ings in Sec. [IV] 


II. SRS RG AS A CONCATENATED BLOCK 
CODE 

The simplest form of real-space renormalization is de¬ 
scribed by a tree tensor network structure, as in Fig. |T| 
From a quantum computing perspective, these tree net¬ 
works are instances of a concatenated block quantum er¬ 
ror correcting code El, and the well-studied properties 
of such codes can provide guidance in choosing an RG 
scheme with an appropriate structure. 

The basic element of an RG scheme of this form is 
an isometry W : Tid' -A Hf b that maps quantum states 
of a single d'-dimensional spin with Hilbert space l~Ld' 
into the Hilbert space of a block of b spins, each with 
Hilbert space Hd of dimension d. The isometry condition 
is H'l W = I d'j with I d' the identity operator on Wd', and 
this requires d' < d b . We consider only schemes with 
d' = d, meaning that a single d-level spin is mapped into 
a block of b identical spins, but it is straightforward to 
extend to the general case d' ^ d. 

In diagrammatic tensor notation, as in Fig. [lj this 
isometry is represented by a tensor with one incoming 


r 1 h Hh Hh H h 


FIG. 1: Real-space RG as a concatenated block code, pre¬ 
sented as a tree tensor network. Here, the block size is 6 = 2, 
and 4 layers of concatenation are shown. 


leg and b outgoing legs: 


W = A (1) 


As an example, a quantum error correcting code that 
encodes a spin (quantum system) into b spins is specified 
by an encoding map given by such an isometry W : TLd —t 



We require that our blocking preserve the on-site sym¬ 
metry, as symmetry considerations have been argued to 
underpin the success of this type of simple real-space RG 
scheme. Tensor network methods for quantum many- 
body systems provides a natural framework to incorpo¬ 
rate symmetries m ■ Let Q be a symmetry group that 
acts on-site through a unitary action U g , g £ Q. That is, 
the symmetry acts on a chain of N spins as JJ ® N , and the 
Hamiltonian H N for N spins satisfies [H N , U® N ] = 0 for 
all g £ Q. Our block encoding then preserves the on-site 
symmetry of the group Q if it intertwines the represen¬ 
tation Uf b on b input spins with a representation V g on 
the output spin, i.e., 

Uf b W = WV g , Vg eg. (2) 

This can be expressed as a tensor relation (shown here 
for a block size of b = 2) as 

I t 

v g G G . (3) 

Ug U g 

Again, we can look to quantum error correcting codes for 
examples. A quantum gate T is said to act transversally 
on a code if the action T® b on all physical spins in the 
block is equivalent to the encoded action of T. Thus, 
any code with transversal action of all gates U gi g € Q 
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provides a blocking that satisfies Eq. ([3| where V is the 
original on-site representation U. 

A real-space RG scheme that has the form of a tree 
tensor network is now entirely determined by identifying 
an isometry W. We now use the structure of the sys¬ 
tem’s Hamiltonian to identify the specific choice of this 
isometry. We seek to separate terms in the total system 
Hamiltonian H into groups H- r n consisting of terms that 
act within a block /, and the remaining terms H lestl such 
that 

H = H Iest + Y. H in ■ (4) 

/£ blocks 

For simplicity, consider the translationally-invariant case 
where all blocks are identical, meaning H( u = H m for all 
I. There is considerable freedom in this decomposition, 
as not all terms with support in a block must necessarily 
be included in H ln . 

Following the approach of Refs. El ED, we place some 
specific constraints on the choice of H ln that will in turn 
yield a suitable RG map. First, note that any choice of 
block Hamiltonian H ln defines a linear projection map 
E 0 : 'Hf b —> Hd 0 , where Hd 0 is the ground space of H ln 
and do its dimension. If do > 1, then the ground state 
space of Hi n is degenerate; from a code perspective, we 
say that such an H ln with a degenerate ground space 
defines a code Hd a C H® b . For our RG scheme, we seek 
to choose our block Hamiltonian H ln with a degenerate 
ground space of dimension do = d, with projection Eq 
onto this degenerate ground space. By defining W = E ' 0 , 
we obtain an isometry as required for our RG map. 

In addition, we require that H ln is invariant under the 
symmetry, [H- ln , U® b ] = 0 for all g £ Q. This condition 
ensures that the ground space 7~Ld 0 of H ln carries a rep¬ 
resentation F of Q. The resulting isometry W = Eq will 
then transform under Q according to Eq. © as required. 
A natural restriction is to the scale-invariant case where 
V = U, the original on-site representation of Q. In this 
case, the degenerate ground space of H ln acting on b spins 
has the dimension of a single spin, and the symmetry 
group Q acts on this degenerate ground space exactly as 
it acts on a single spin, i.e., transversally. 

One round of blocking for our RG map renormalizes 
the Hamiltonian as H —> FT, where 

H’ = (<SjW / t ).ff(®jWj) . (5) 

As by construction we have that Wj H J m Wj = cld for 
Id the identity operator on the renormalized spin from 
block I and c a constant, the renormalized Hamiltonian 
is determined by 

H’= {®iW\)H iest (® I W I ), (6) 

up to an irrelevant additive constant, which is now an 
operator acting on the renormalized spins. 

A real-space RG scheme such as this has been applied 
to the transverse field Ising model by Miyazaki et al. [ 8 ] , 


based on the Z 2 spin-flip symmetry of the Ising model, 
where the resulting 2-to-l RG produces the exact expres¬ 
sion for the critical exponent v that describes the scaling 
of the Ising model correlation length. This method is also 
surprisingly accurate when applied to 2D quantum Ising 
models on a variety of lattices ElE]- 

III. REAL-SPACE RENORMALIZATION OF 
THE QUANTUM ASHKIN-TELLER MODEL 

In this section, we will apply the SRS RG method de¬ 
scribed above to the Ashkin-Teller model. Note that the 
critical behavior of this model has previously been in¬ 
vestigated using real-space RG techniques in Ref. ESI, 
wherein both a standard blocking method and a related 
(but distinct) symmetry-respecting method were used. 

A. The quantum Ashkin-Teller model 

The quantum Ashkin-Teller (AT) model [TOi [TTj is de¬ 
fined by the Hamiltonian 

Hat = — J YXZ^Z^, + Zj^Zj+ 1,2 

3 

+XZj t iZj + i t iZj t 2Zj + i > 2) 

-hY(X Jtl + X j:2 + XX jtl X Jt2 ) , (7) 

3 

with J, h, and A real-valued coupling coefficients. This 
model can be viewed as a pair of transverse-field Ising 
chains coupled by two- and four- spin terms. In partic¬ 
ular, when A = 0, we have a decoupled pair of quantum 
transverse-field Ising chains. 

The AT model possesses an on-site D 4 symmetry, de¬ 
fined as follows. Define the operators 

Si = 0 X jt i , 5 2 = 0 X jt 2 . (8) 

3 3 

These operators commute with the Hamiltonian, and 
generate a Z 2 x Z 2 symmetry. The Hamiltonian is also 
invariant under the symmetry that swaps the two chains, 
i.e., Xj t 1 -f-> Xj 2 and Zj \ -f-)- Zj <2 . Together, these sym¬ 
metries form the symmetry group D 4 . (We note that this 
model also possesses a non-local symmetry, given by ap¬ 
plying the self-duality map of the Ising model to each of 
the two chains. This self-duality ensures that, for a fixed 
value of A, if there is a single phase transition it must be 
at h = J.) 

The model is critical along the line /3 = h/ J = 1 
and — 1 /a/ 2 < A < 1, and the critical indices vary con¬ 
tinuously along this line. Based on finite size simula¬ 
tion EU, CFT arguments [18], and real-space RG based 
on MER A [13], the critical line of the AT model has been 
identified as the orbifold boson with radius Ro given by 

R b = |[arccos(-A)] _1 , (9) 
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for A € (—1 /a/ 2, 1] - In addition, the AT model for 
/3 = 1 for any A can be transformed into the six-vertex 
model [12], and so is exactly solvable using Bethe Ansatz. 
These relations along the critical line to the six-vertex 
model and the orbifold boson CFT allows us to obtain 
analytic predictions for the critical behavior of this model 
along the line f3 = i mm, which is governed entirely 
by the radius Rq of the orbifold boson. 

In particular, consider the correlation length £ defined 
as for the Ising model — (Z it i} 2 ) ~ e -0'-*)/S j 

where we have arbitrarily chosen the first Ising chain but 
symmetry ensures that it is the same for both chains. 
With /3 = h/J as a tuning parameter and /3o = 1 at 
the critical point, the correlation length near the critical 
point is governed by £ = (/3 — /?o) _1/ *, with v * the corre¬ 
lation length critical exponent. This critical exponent is* 
is expected to vary continuously along the line of criti¬ 
cality. From Ref. [12], we have a theoretically predicted 
value of j/* given in terms of the orbifold boson radius 
Ro as v * = 1/(2 — R 2 q ) and therefore 

V * 2 — |[arccos(—A)] -1 

Based on the success of predicting this critical expo¬ 
nent for the quantum transverse-field Ising model (where 
z/*(A = 0) = 1) using the SRS RG, we now aim to cal¬ 
culate a prediction for this exponent in the quantum AT 
model. 


B. Symmetries and block Hamiltonians 


a 4-dimensional degenerate ground state on which the 
D 4 symmetry acts transversally. These properties are 
seen most easily from the perspective of the Hamiltonian 
code defined by this block Hamiltonian’s ground space. 
Defining the logical operators for a 2-spin phase-flip re¬ 
dundancy code 

A’i = Xj ^Xj + i^ , Zi = Zj + 11 , (12) 

X 2 = Xj 2 Xj + 12 , Z 2 = Zj - 1 - 1,2 j (13) 

we see that these operators commute with H ul , acting as 
a pair of renormalized Pauli operators on the degenerate 
ground space. 


C. Solving the block Hamiltonian 


To find an explicit expression for the isometry W de¬ 
fined by an embedding of the encoded site into the ground 
space, we need to solve H ln . For convenience, we change 
to a new set of operators 

A) = -Xj'.i; Zi = Zj iZj + i i , (14) 

*2 = X jt 2 , Z 2 = Z ji2 z j+li2 . (15) 


These operators commute with the logical operators of 
Eq. (41), and we can express H ln entirely in terms of 
these operators as 


i„ = ~Z X -Z 2 - pX! - pX 2 

— \/3X\X 2 — \Z\Z 2 , (16) 


Here we describe an SRS RG map for the Ashkin-Teller 
model as a concatenated block code that respects the D 4 
symmetry of the model. This on-site -D 4 symmetry acts 
on a pair of spins (labelled by (j, 1 ) and (j, 2 ) in the 
Hamiltonian of Eq. 0). that is, a ‘site’ consists of a pair 
of spins and is labelled by j. 

The simplest blocking, which we consider first, is a 2- 
to-1 blocking, i.e., that maps 4 spins (two sites) to 2 spins 
(one site). We therefore seek an isometry W : H 4 —> 
Tif 2 that encodes a 4-dimensional quantum spin into a 
pair of 4-dimensional spins. To identify an appropriate 
isometry, we select terms in the Hamiltonian ([T]) to define 
a block Hamiltonian H in on a pair of sites, possessing a 
4-dimensional degenerate ground state on which the D 4 
symmetry acts transversally. 

As a starting point based on the successful solution of 
this form for the quantum Ising model with its Z 2 sym¬ 
metry as in Refs. El, one choice of a block Hamiltonian 
satisfying these conditions is 

Hln = —J(Zj,lZj+i t l + Zj 2 Zj + i 2 

+ \Zj^iZj + i.iZj^ 2 Zj + i t2 ) 

- h(X jt 4 + X jt 2 + A.v, ,.V ; . 2 ) , (11) 

expressed here as acting on sites j and j + 1 within a 
block (e.g., by choosing j odd). This Hamiltonian has 


where we recall that /3 = h/J. We note while H ln is 4-fold 
degenerate on the 4 spins in the block, the H amil tonian 
Hi n is nondegenerate in terms of the operators (15). This 
is because H ln is supported on a single pair of effective 
spins, as follows. It is most natural to view the space 7if 2 
of the two sites as possessing a virtual tensor product 
structure ^ 407 / 4 , where R 4 corresponds to the support 
of the operators (15) and H 4 the support of the logical 
operators (41). Thus, the block Hamiltonian takes the 
form H ln = H m < 8 > I 4 -, where having H ln nondegenerate 
ensures that H in is 4-folcl degenerate. 

This Hamiltonian can be solved analytically. It has a 
particularly simple solution along the line ;3 = 1, where 
the ground-state energy is given by 


E g (P = 1, A)/J = —A — \/A 2 + 8 , (17) 

and the ground state \gp=i,\) can be expressed as 

\ 9 p-i,\) = cos( 0 / 2 )| <j>)\ 4 >) + sin( 0 / 2 )|</- L )|</r L ). (18) 


Here, we have defined the state | <f>) as the +1 eigenstate 
of the operator (X + Z)/y/2 , | (j)- 1 ) is the orthogonal spin 
state with eigenvalue — 1 , and tan 0 = A/(2\/2)- 

The general solution for the eigenvalues and eigen¬ 
states of the block Hamiltonian H ln of Eq. (JTgJ) are closed- 
form cubic expressions (not shown). The ground state 
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\gp,x) is nondegenerate for A £ (—1,1], and there are no 
crossings. 

The isometry W then takes the simple form 

W = \gp,x) ® h , (19) 


in terms of the virtual tensor product structure H 4 ® H 4 
given by the operators (15) and (41). 


D. Renormalizing the Hamiltonian 


1. Identifying the critical line 


From Eqs. ( 23p4 | and our solution | gp,x) , we find a line 
of fixed points for the h and J coefficients along the line 
(3 = h/ J = 1. That is, from the flow of these coefficients 
we recover the known critical line. 

We note, however, that the coefficients A x,z are not 
fixed points along this entire line for the SRS RG map. 
We can use our analytical solution of Eq. (181 along the 
critical line /? = 1 to evaluate 


Applying one round of blocking for our RG map results 
in a renormalized Hamiltonian given by 


H' = (® J W 7 t )H r (0/W / ) , 


( 20 ) 


which, in terms of our new operators and the solution 
I^a) for the ground space of H in , takes the form 


H' = - 




7 ^ 1 , 7+1 


+ {gp,x\Z2\~gp,x)z 2J z w 
+ ^{gp,\\ZiZ 2 \gp'\)ZijZ 1 j + iZ2jZ2j+i 
-h'^2^[gp,\\Xi\gp'\)Xij + {gp,x\X 2 \gp i x)X 2 j 


+ M9p,x\XiX2\gp,\)Xi i iX2,i 


( 21 ) 


where / denotes the I th site in the renormalised chain (a 
block of two sites in the original chain). 

Because the Hamiltonian H- ln has a symmetry swap¬ 
ping 1 -f-)- 2, we have identities (gp,x\Zi\gp t x) = 
(gp,\ | Z 2 \gp,x) : etc. Therefore, the renormalized Hamil¬ 
tonian can be expressed as 


■s 

II 

,iZij + i — J 1 Z 2 jZ 2 j+i — h'Xij — 

h'X 2J 

- X' x h'X 

ijX 2 j — X' z J' ZijZij + iZ 2 jZ 2 j + i 

(22) 

where the coefficients are given by 


J'/J = 

{gp,\\Zi\gp t x ), 

(23) 

h'/h = 

(gp,x\Xi\gp,x), 

(24) 

^xt'Xx = 

(h/h')(gp t x\XiX 2 \gp^x) 


= 

(gp,x\XiX 2 \gp t x)/(gp,x A'i .9/3 ,a) , 

(25) 

A' z /A z = 

(J/J')( 9 p,\\ZiZ 2 \gp t x) 


= 

(gp,x\ZiZ 2 \gp,x)/(gp,x\Zi\gp : x) ■ 

(26) 


We note that the two terms with A in the Hamiltonian do 
not rescale the same way, and so we have defined Xxh as 
the coefficient of the XX term and A zJ as the coefficient 
of the ZZZZ term. The symmetry of the AT model that 
relates these two terms is the self-duality symmetry of 
the Ising model, which is a nonlocal symmetry that is 
not enforced in the SRS RG map. 


{gp,x\Zi\gp,\) =-j==^ (27) 

{gp,\\X\\gp : \) =-j=== (28) 

(~gp,x\Z 1 Z 2 \~gp, x ) = ±(l + -^=) (29) 

{gp,x\XiX 2 \gp,x) = - (l + y=j=|) • (30) 


We then find that the coefficients rescale along the (3 = 1 
line according to 


h' = 
J' = 
A'= 


. = h 

VA2T8 

2 J 

vra 

^(a + a/wTs). 


(31) 

(32) 

(33) 


Note that the SRS RG map has two fixed points, A* = 0 
and A* = 1. Note these fixed points are the same as those 
found in the related real-space RG scheme of Ref. [16]. 

When the SRS RG map is evaluated along the critical 
line /3 = 1, the two parameters Xx and A z flow in the 
same way, i.e., the expressions (251 and (26) are identical 
at /? = 1. (This is not true for general pi] As a result, 
when considering the critical line we can dispense with 
the need for two different coefficients and return to using 
the single coefficient A to describe critical behavior. 


2. Calculating critical exponents 

We now use the SRS RG map to determine the be¬ 
havior along along the critical line j3 = 1. Although the 
SRS RG map results in a flow along this critical line to 
two fixed points (A* = 0 and A* = 1), we show that we 
can use a linearised solution around any point along the 
model’s critical line to estimate the correlation length 
critical exponent of the AT model, which turns out to be 
independent of the SRS RG map’s flow along the critical 
line. 

About any point along the critical line (/3 = 1, A £ 
(— 1 /a/ 2, 1]), the behavior is governed by the Jacobian of 
a two-parameter transformation (/3, A) —> (/3',X'). The 
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closed-form solution for the matrix elements of this Ja¬ 
cobian using Eqs. (23p6) and \gp,\) is 


df¥_ 

d/3 

av 

d/3 

d/T 

dA 

av 

dA 


3=1 


= 2+^(-A+VA 2 + 8) 


3=1 


3=1 


3=1 


= 0 , 

— g" — A“)(—A + \/A 2 + 8), 

_ 2 + i\(\ + y/)?+8) 

VF+8 


(34) 

(35) 

(36) 

(37) 


The eigenvalues of the Jacobian, which can be expressed 
in the form b Vl ’ 2 where b is the block decimation size 
(5 = 2 in this example, being a mapping of 4 spins to 
2), and v\^ 2 are functions of A € (— 1 /a/ 2, 1]- As the Ja¬ 
cobian is a triangular matrix, the eigenvalues are simply 
the diagonal entries, i.e., 


2 /i (A) — log 2 2 + — (—A + \J\ 2 + 8) 


2/2 (A) = log 2 


2 + ±A(A + v'A r +8) 1 


y/W+8 


(38) 

(39) 



FIG. 2: Comparison of the critical index v* of the energy 
operator (black) with the critical exponent v obtained from 
our real-space RG map (red) as functions of the Hamiltonian 
parameter A, using the smallest blocking size b = 2. Inset: 
Relative error \v — 14 I /14 of the predicted critical exponent 
v obtained from our real-space RG map to the critical index 


The first of these eigenvalues y\ associated with 
d/3'/d/3 gives the estimate of the correlation le ngth c riti- 

the 


III A 


cal exponent as v = . As discussed in Sec. 

value of this critical exponent is predicted to be v * as 
given by Eq. ( |l0| ). Comparing i/(A) and 14 (A) in Fig. [2j 
we see broad agreement, including exact agreement at 
A = 0, where i/(0) = 14 (d) = 1 but decreasing accuracy 


towards the limits A = — 1/V2 and A = 1. For com¬ 
parison, at A = 1 we have 14 ( 1 ) = 2/3 compared with 
i/(l) = (log 2 5 — l ) -1 ~ 0.75. (Note that, at A = 1, the 
AT model is equivalent to the 4-state Potts model.) At 
A = —l/y/2, we have 14 (—l/\/2) —► 00 compared with 
v{-\/y/2) - 2.25. 

The second eigenvalue 2/2 is associated with the flow 
along the critical line that is predicted by the SRS RG 
map, with fixed points A = 0 (with y 2 = —1/2 at this 
point) and A = 1 (with y 2 = { 2 — log 2 3)). We contrast 
this flow to predicted behavior of the AT model, where 
the corresponding scaling field is marginal. 


E. Alternative blockings 

Our analysis so far has restricted to blocking of 4 spins 
to 2 (i.e., 5 = 2). In this section, we explore alternate 
blocking choices for the SRS RG of the quantum AT 
model. 

For an alternative blocking choice b , we seek an isome¬ 
try Wb : Ha —>• Hf b that encodes a 4-dimensional quan¬ 
tum spin into b such spins. Following Sec. EH a natural 
generalization of the 5 = 2 case is to use the phase-flip 


redundancy code with logical operators 


Xi — Xj'iXj+i'i ■ ■ ■ Xj + b- 1 , 1 , Z\ — Zj+k 0 , 1 5 (40) 

X 2 = -^0, 2 -^+ 1,2 • • • Xj+b- 1,2 , Z 2 = Zj+k Q 2 1 (41) 

where ko £ {0,..., 5 — 1} is an integer specifying a pre¬ 
ferred site in the block. These operators commute with 
the following choice of block Hamiltonian 


b—l 

Hin = -j£(Z,. +fe , iZj+k+1,1 + Z j+k , 2 z j+k+1 , 2 

k =0 

+ ^Zj+k,lZj+k+l,lZj+k,2Zj+k+l,2) 

— h [Xj+k+ + Xj+k,2 + A Xj+k,iXj+k,2) ■ (42) 

Note that this block Hamiltonian includes all terms in the 
quantum AT Hamiltonian with support entirely within 
the block, except for the X- and A'V-type terms on the 
preferred site ko. This choice of block Hamiltonian has a 
4-fold degenerate ground space, and therefore defines an 
isometry W with the desired properties of our real-space 
RG map. 

We consider two choices of preferred site, closely fol¬ 
lowing Ref. PJ: ko = 0, corresponding to the preferred 
site at the edge of a block, and k 0 = (5 — l)/2 for odd 5, 
corresponding to the middle site of a block. For 5 > 2, 
we solve H m and determine the isometry W numerically, 
for a range of block sizes as illustrated in Fig. [3] 

The predicted critical exponent v is found to be largely 
independent of the size of blocking 5, as shown in Fig. [4] 
All choices of blocking recover the exact result v* = 1 
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b = 2, ko = 0: 
b — 3, ko = 0: 
b — 4, ko = 0: 
b = 5, fco = 0: 
b = 3, fco = 1: 
b = 5, k 0 = 2: 




FIG. 3: Alternative choices of blocking. Within each block 
consisting of b sites, one site (pair of spins, such as shown 
encircled by dashed line), shown in red, is a preferred site for 
the encoding map. 


at A = 0, as in Ref. [S], For A > 0, all blockings with 
preferred site ko = 0 at the edge give nearly indistin¬ 
guishable predictions, where for A < 0 the accuracy in¬ 
creases with increased blocking size but only marginally, 
with negligible improvement observed beyond 6 = 4. For 
6 = 3, we show a comparison of the choice of preferred 
spin k 0 = 0 (edge of the block) and k 0 = 1 (center of 
block), shown in Fig. [5j Again, both cases recover 14 = 1 
at A = 0 exactly. The choice ko = 0 gives more accurate 
results for A > 0. For A < 0, we see the largest qual¬ 
itative differences of all blocking choices but still have 
broad agreement. For 6 = 5 we have performed a similar 
comparison of ko = 0 and ko = 2 (not shown), where we 
find similar but less pronounced differences. 


IV. CONCLUSIONS 

We have presented a generalization of the simple real- 
space RG scheme proposed in Ref. jS] and further ex¬ 
plored in Ref. [9], based on a blocking isometry that re¬ 
spects the symmetry of the model. We have used this 
SRS RG to investigate the critical behavior of the quan¬ 
tum Ashkin-Teller model, which possesses a line of crit¬ 
icality along which the critical exponents vary continu¬ 
ously, and which reduces to decoupled Ising models at 
one point on this line. The SRS RG map was used to 
extract estimates of the correlation length critical expo¬ 
nent v as a function on this line of criticality for various 
blocking schemes, including a closed-form expression for 
the simplest choice of blocking. 

The SRS RG identified the line of criticality in this 
model, and the resulting estimate for the correlation 
length critical exponent v varied continuously along the 
critical line. Compared with the exact quantity j/*, it 


FIG. 4: Critical exponent v obtained from our real-space RG 
map using various blocking sizes b, all with the preferred site 
ko = 0 on the edge of the block, as illustrated in Fig. [3] For 
comparison, the exact value 14 is shown (black line). Inset: 
Relative error \v — v*\/v* of the predicted critical exponent v 
obtained for each choice. 



FIG. 5: Critical exponent v obtained from our real-space RG 
map using a blocking size 6 = 3 with a symmetric (orange) 
vs asymmetric (blue) choice of internal Hamiltonian terms as 
illustrated in Fig. [3] For comparison, the exact value is 
shown (black line). Inset: Relative error |v — of the 

predicted critical exponent v obtained for each choice. 


gave a precise value at the Ising decoupling point and re¬ 
produced the broad features of its functional form along 
this line as a simple analytical expression with the small¬ 
est block size. In a range of A around the decoupled 
Ising point A = 0, the accuracy of the SRS RG is good, 
with a relative error in the prediction of < 10 % in the 
range A £ [—0.4, 0.9]. We compare this accuracy to that 
obtained by the sophisticated and numerically-intensive 


































































MERA simulations performed in Ref. Ha. where the rel¬ 
ative error for this critical exponent was ~ 5% over this 
same range. Exact diagonalization (also in Ref. ]T3]) gave 
a comparable error to MERA. 

The estimate v became an increasingly inaccurate pre¬ 
dictor of v* towards the ends of the critical line. We 
note that, in a previous study of the critical behavior 
of the quantum AT model using MERA [13] (which is 
also a form of real-space RG), the predicted critical in¬ 
dex associated with this exponent also deviated upwards 
from the exact value towards the endpoint A = 1 of this 
critical line, in a way that is very similar to the results 
shown here. Exact diagonalization calculations produced 
a similar behavior. 

As noted previously, this simple approach to real-space 
RG can accurately reproduce some of the critical expo¬ 
nents of a model, but not necessarily all of them [9]. For 
the Ising model, the scheme of Ref. 0 accurately pre¬ 
dicts the correlation length exponent, but not the mag¬ 
netic exponent. In the quantum AT model, there is a 
rich critical behavior associated with the scaling dimen¬ 
sions of the orbifold boson CFT, included twisted scaling 
dimensions that do not appear in the standard c = 1 bo¬ 
son CFT. Although such critical indices are not observed 
here (the lowest such index is 1/8, independent of A), the 
success in MERA for identifying these additional critical 
indices |T3] gives hope for generalizations of simple real- 
space RG structures such as the one presented here to 
predict further aspects of the critical behavior. 

The SRS RG method can yield a simple, analytical 
RG map based only on the symmetry properties of the 


model and yet can also qualitatively reproduce aspects 
of the critical behavior. As such, it may be useful as a 
simple approximation method even if its accuracy is not 
as competitive as existing numerical tools. In particular, 
it may be useful is assessing the role of different symme¬ 
tries that may be applied to a class of models [15]. We 
note that although the quantum AT model studied here 
is characterised by an Abelian Z 2 x Z 2 symmetry, the 
SRS RG applies quite directly to continuous symmetries 
such as SU(2). While the intertwining requirement of 
our RG map given in Eq. ®> can give a much broader 
range of solutions for such continuous symmetries, we 
note that for characterising zero-temperature phases for 
spin chains with continuous symmetries such as SU(/c) 
and SO(2fc + 1), it is sufficient to characterise the rep¬ 
resentations of an Abelian subgroup (in these cases, the 
subgroups (Zk x Zk) and respectively) [20]. 
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